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Abstract 

Classical scalar fields have been proposed as possible candidates for the dark matter component 
of the universe. Given the fact that super-massive black holes seem to exist at the center of most 



OO ^ 

. galaxies, in order to be a viable candidate for the dark matter halo a scalar field configuration should 



be stable in the presence of a central black hole, or at least be able to survive for cosmological time- 
scales. In the present work we consider a scalar field as a test field on a Schwarzschild background, 
and study under which conditions one can obtain long-lived configurations. We present a detailed 
study of the Klein-Gordon equation in the Schwarzschild spacetime, both from an analytical and 
numerical point of view, and show that indeed there exist quasi-stationary solutions that can 
remain surrounding a black hole for large time-scales. 
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I. INTRODUCTION 



Scalar fields are a common theme in modern cosmology. They play a central role in 
inflation and they have been frequently used to describe dark energy in place of the 
cosmological constant {2, 3|. Moreover, even if the usual description of dark matter is done 
in terms of weakly interacting massive particles (WIMPs), given the persistent uncertainty 
about the real nature of that mysterious component, one can not a priori reject one candidate 
in favor of another. In this context, ultra- light scalar flelds have also been invoked as possible 
alternative dark matter candidates j4-21|. 

In modern physics quantum flelds are considered to be the fundamental constituents of 
nature. For the particular case of a boson fleld which has not been in thermal equilibrium 
with the other flelds of the standard model of particle physics, a classical fleld description 
(i.e. a coherent excitation) seems in principle as natural as that given by the WIMP hy- 
pothesis. Indeed, a description in terms of classical scalar flelds has been shown to describe 
the large scale cosmological scenario (to linear order) as consistently as cold dark matter 
(CDM) [ill [igI 22, 23|. In a more local context, massive complex scalar flelds can form 



spherically symmetric, asymptotically flat, regular, self-gravitating conflgurations known as 



25| (see also 26|, l27| and references therein). For the case of a scalar fleld 



"boson stars" [24, 

with a mass of the order of 10~^^ - 10~^^eV, these conflgurations can have the mass and the 



size expected for a typical galactic halo J5 
though such conflgurations are stable (291 . 
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However, it has been recognized that even 



30| . they do not reproduce the observed velocity 



proflles (though this problem can be signiflcantly alleviated with the use of multi-state con- 



flgurations, see e.g. 



31 



-|33| ) . Finally, ultra-light scalar fleld dark matter (SFDM) could have 



some advantages over the standard Lambda cold dark matter (ACDM) scenario, narnely that 



the SFDM model is not expected to produce an over-density of satellite galaxies 



11| . nor 



does it result in the generation of the cuspy halos typically obtained in the high-resolution 



2l|. 



N-body simulations performed within the standard ACDM model 

Despite all this, there has been some controversy regarding the possibility that dark 
matter halos could be described by coherent scalar excitations. The central issue being the 
fact that most galaxies seem to host a super-massive black hole at their centers, so that a 
scalar fleld halo should be able to survive surrounding such a black hole for cosmological time- 
scales. Regarding that question, it has been known for a long time that static, spherically 



symmetric, asymptotically flat scalar field configurations with a Schwarzschild black hole 



at their center can not exist in nature 
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"scalar hair" 



34j . i.e. a Schwarzschild black hole can not have 



However, these "no-hair" theorems have nothing to say about a specific 



time-scale, and it could in principle be possible that scalar field configurations could survive 
the birth of a super-massive black hole and remain in a quasi-stationary state for very long 
times. 

The question of how a scalar field is accreted into a Schwarzschild background has at- 
tracted the attention of many researchers, motivated in some cases by the fact that dark 



energy 



els 



44 



can be modeled by a scalar field 



37H43| . or because of the scalar dark matter mod- 



45] . A definite picture of the accretion of a scalar field onto a black hole is not clear 



yet and there is in fact some controversy around it. For instance, in [4^ it is argued that the 
accretion rate will be very slow, while others argue that this accretion will be fast 45|. Very 
long lasting configurations have also been found using "artificial" scalar field potentials 46 |. 

Here we will restrict ourselves to the case of a scalar field as a test field in the background 
spacetime of a Schwarzschild black hole. It is for this reason that we will talk about "scalar 
field clouds" rather than "galactic halos", in order to maintain a clear perspective of the 
real problem we are studying. The dark matter content in the central region of a galactic 
halo is very diluted, so it is to be expected that this approximation could shed some light 
on the nature of the dark matter problem. However, we would like to emphasize that this 
is only a first step on the study of long-lived scalar field configurations surrounding a black 
hole, and the self-gravitating case will be considered in future works. 

This paper can be separated in two distinct parts plus the conclusions. In a first part 
(Section [IT]), we introduce the problem and present some important analytical results. We 
find that it is indeed not possible to obtain a stationary, localized, physical scalar field 
configuration surrounding a Schwarzschild black hole: Solutions that are well behaved at 
spatial infinity are not regular at the horizon; they oscillate with a divergent frequency 
and these oscillations produce an infinite integrated energy close to the horizon. Thus, in 
particular, such solutions are not compatible with the test field limit approximation. This is 
of course not surprising if we remember the no-hair theorems. However, these non-physical 
solutions can be quite useful in order to construct appropriate initial data for a special set 
^ Here we will be working within the domain of general relativity. However, "hairy" black hole eonfigurations 
are indeed possible in the context of scalar-tensor theories of gravity 36| . 
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of physical quasi-stationary (resonant) scalar field configurations with very long lifetimes. 
In a second part (Section IIIII) . we perform numerical evolutions in order to determine the 
lifetime of such quasi-stationary resonant configurations. Finally, we discuss our findings in 
Section dVl 

Throughout the paper we will use geometric units such that G = c = 1, and we will use 
the conventions of Misner, Thorne and Wheeler (MTW) In particular, the signature 
of the spacetime metric will be taken to be (— , +, +, +). 



II. EIGENMODE ANALYSIS AND RESONANCES 



We will be looking for stationary scalar field configurations surrounding a Schwarzschild 
black hole. Here we will restrict ourselves to the case of a canonical, massive, non-self- 
interacting, minimally coupled scalar field 0. The equation of motion for that field is given 
by the Klein-Gordon equation 

(□-/i')0 = O, (1) 

with the d'Alambertian operator defined as □ := (1/ y/—g) dfj^{y/—gg^^dy). With our con- 
ventions is dimensionless, while /i has dimensions of (length)"^. The associated quantum 
mechanical "mass" of the scalar field is given by fl/i, but all our considerations in this 
work will be purely classical. (For the SFDM models usually considered in the literature 
one usually has ~ 10~^^ - 10~^^eV in physical units, which implies /i of the order of 
1 - 10-2(pc)-i.) 

We will also be assuming that the energy associated with the scalar field configuration is 
very small compared to the mass of the black hole, so that the gravitational back-reaction 
associated to the scalar field distribution can be disregarded. Under this assumption, and 
working in Schwarzschild coordinates, we can write the spacetime metric as 

ds^ = -N{r)dt^ + + r^dfi^ ^ jY(^) ^ _ ^Mjr , (2) 

with M the mass of the black hole and d^ := dd'^ -|- ^m^ddi^'^ the standard solid angle 
element. This coordinate system is suitable to describe the physics in the exterior region 
r e (2M, oo). 

To look for solutions of the Klein-Gordon equation above, we will start by considering a 
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decomposition into spherical harmonics: 



(3) 



£,m 



where the 1/r factor has been introduced for convenience, and the parameters i and m take 
the usual values: i G {0, 1, 2, . . .}, —£< m < i. Substituting the above equation into ([1]) we 
obtain the following family of reduced equations 

1 



im 







where we have defined 



+ 1) 2M 



^2 ■ ^3 • 



+ 



(4) 



(5) 



Notice that is an everywhere positive and monotonically decreasing function of r, and 
also that equation (jl]) does not depend on m. The first term in ([5]) is the usual centrifugal 
term, the second one is a curvature correction, and the third one is inherited from the 
corresponding term in the Klein-Gordon equation, Eq. ([1]). 



In the test field limit all configurations have the conserved energy E = ^ Eim, with 

£,171 



£m 



pE{r)dr 



2M 



and where 



pE[r} 



N{r) 



£m 



dt 



+ N{r) 



di'£„ 



dr 



(6a) 



(6b) 



is the energy density integrated over the sphere as a function of areal radius. In what follows 
we will refer to this integrated energy density simply as the "radial energy density" . 

To obtain the last expression we have used the fact that it is always possible to associate 
a conserved charge Qa = ka^Tf^^n^y/^ dT, to every Killing vector field ka- Here is 
a time-like, future-directed, normalized 4-vector orthogonal to the 3-dimensional space-like 
hypersurface S, and y/^dT, is the volume element on S. The subindex a labels the different 
Killing fields associated with the background metric, and T^^ is the energy-momentum tensor 

= (— ml i-^- iiot all the modes iperiiit,^) in (jSj are 



^ If the scalar field was real we would have ipi.-' 
independent of each other. Apart from this, all the results obtained below apply for the real as well as the 



complex cases. The differences between the real and complex scalar field will appear when considering the 
energy-momentum tensor, but this is not relevant for the test field limit approximation we are considering 
here. 
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for the matter fields. Here we will be only interested in those configurations with a finite 
value for the integrals above. Under this assumption and the additional requirement of S 
being a Cauchy surface, it follows that the charges Qa do not depend on the particular choice 
of S; in particular the charges are preserved under time evolution. For the case = dt 
and S = {t = const.}, where {t,r,6,(f) are the usual Schwarzschild coordinates, we arrive 
at expressions ([6]). The hypersurfaces S are Cauchy surfaces for the Schwarzschild exterior 
spacetime, and as a consequence the energy ([6]) is conserved. However, in the numerical 
simulations of Section IIIII below, we will be integrating the energy density over a different 
foliation which extends from the future horizon to some finite, large radius, and the resulting 
energy will, in fact, decay due to energy loses through the horizon and the outer boundary. 



A. Reduction to a time-independent Schrodinger-like problem: mode solutions 

In order to look for the stationary solutions of Eq. (jlj) we make a further decomposition 
of the functions ipirn{t,r) into oscillating modes of the form: 



(7) 



with ugm a real frequency, and Uim{r) a complex function of r in the interval (2M, oo). In 
order to simplify the notation from now on we will omit the labels £ and m from uem and 

Introducing the ansatz ([7]) into equation (jl]) one finds the following eigenvalue problem 



-iV(r)|- (iV(r)|: ]+N{r)Ueifi,M-r) 



u{r) = uj'^u{r) , 2M <r <oo . (8) 



After defining the Regge- Wheeler tortoise coordinate r* := r + 2Mln(r/2M — 1), this equa- 
tion can be rewritten as the following time-independent Schrodinger-like equation: 



dr 



*2 



uir*) = Lo'^uir*] 



-oo < r < oo. 



(9) 



with the "effective potential" Vcs{r*) defined as 



Ves{r*) := N{r)Ue{fi,M;r) , r = r{r*). 



(10) 



Notice that in the usual non-relativistic quantum mechanical case one would have u instead 
of cj^ in Eq. ([9]). This is because the Schrodinger equation is first order in time, whereas the 
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FIG. 1. The effective potential M'^V^s for i = 1, and different values of the parameter M/i. The 
shaded regions represent the "resonance bands" where the potential has a local minimum (see 
Section He] below). Notice that for Mfi = and Af/i > 0.466 the local minimum does not appear. 
For the case Mfi = we recover the usual Regge- Wheeler potential. The case with Mfi = 0.3 will 
be considered later on in our future examples and it appears as a solid line to emphasize it. 



Klein-Gordon equation is second order. Also, note that equation ([9]) in fact only depends 
on the parameters i, Mfi and Mu. Since that equation is invariant under u — —u we can 
restrict ourselves to the case a; > 0. 

The typical form of the effective potential can be seen in Fig. [H Contrary to what happens 
in the usual Regge- Wheeler equation, the asymptotic value of V^s for large radii is not zero 
but /i^. This fact will turn out to be crucial in the study of the behavior of the solutions to 
equation ([9]). 

In the following, we will analyze the spectrum of the Schrodinger operator given above 
in more detail. The operator can be separated into the point and the continuous spectrum. 
The point spectrum consists of the eigenvalues of Eq. ([9]) for which the square norm 

+ Ui{fx,M;r)\u\^^ dr (11) 
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uf ■= [ ( iV(r) 
J2M \ 



is finite. Such solutions, when they exist, describe bound states in the quantum mechanical 
picture. In the context of the Klein-Gordon equation they will give rise to stationary modes 
ip{t,r) = e*'^*M(r) with finite energy. However, as we will show below, such solutions do not 
exist in our case. 

This leaves us with the continuous spectrum consisting on the values of for which the 
time-independent Schrodinger equation ([9]) admits solutions which do not grow faster than 
a polynomial as r* — )■ ±00. These solutions describe scattering states and will be analyzed 
in more detail below. 

At this point it is important to mention that exact solutions to equation are in fact 



known and can be expressed in terms of the so-called "Heun functions" (see e.g. 48|, |49| ) . 
However, we will not use such exact solutions here and will rather concentrate on the general 
behavior of the solutions. 

We will start with an analysis of the asymptotic behavior of the solutions to Eq. both 
close to the horizon and close to infinity. 



B. Asymptotic behavior and spectrum 

Regarding the behavior of the effective potential notice first that, as already mentioned 
above, in the limit r* — )■ 00 it approaches the constant value /x^, so that Eq. Q reduces to 

^'"^ +k\ = 0, e := (12) 



*2 



dr 

The general solution to this equation is a linear combination of the functions u = e^"^^^ , with 
k := >y u"^ — fi"^. Therefore, as long as < < /i^ the parameter k will be purely imaginary 
and we will have one exponentially decaying and one exponentially growing solution, while 
for > /i^ the parameter k will be real and both solutions will be oscillatory. 

For r* — 7- —00, on the other hand, the effective potential Vcs decays exponentially fast to 
zero since N{r) = 2Mr^^{r /2M — 1) = 2Mr~^exp((r* — r)/2M). Therefore, the solutions 
are oscillatory for all u"^ > (notice that we are assuming u to be real so this is always the 
case), and they can be expressed as a linear combination of the functions u = e"*^*"^*. 

Consider now the interval < u"^ < fi^. For u"^ to belong to the spectrum, the right 
asymptotic solution must be proportional to e"''"'!''*, and the left asymptotic solution may 
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be a linear combination of e^*'^''* of the form 

u oc e"^'^"* + Re'"^^* , (13) 

with R a reflection coefficient. However, since the right asymptotic solution and the co- 
efficients in Eq. ([9]) are real, so must be the left asymptotic solution, which implies that 
|-R| = 1 such that u oc cos(a;r* — 5) for some phase factor 5. In particular, the solution is not 
normalizable with respect to the square norm given in fllip . and so it does not give rise to 
a bound state. Therefore, the spectrum for < < /i^ is continuous and non-degenerate. 
From the point of view of scattering theory, = 1 simply means that a plane wave which 
travels from r* — )■ — oo to the right with frequency < < /i^ is completely reflected 
from the potential since it does not have enough "energy" to reach r* — +oo. For the 
Klein-Gordon equation, this means that the corresponding stationary mode solution satis- 
fles ip{t,r) oc e*'^(*~^'*) + i?e*'^^*"'"^'*-* for r* — )■ — oo. Since i? 7^ 0, the solution is always singular 
at both the future and past horizons. 

Finally, for u"^ > fi"^ the asymptotic solutions at r* — ±00 are both oscillatory, and the 
spectrum is continuous and degenerate since there are two linearly independent solutions 
which are bounded at r* — > ±00. 

Notice that the mode solutions with < < /i^ do not have a flnite value for the 
energy, expression (|6]), but they are still "localized" around the black hole since they decay 
exponentially at large r*. This is not the case for the mode solutions with co^ > n"^, which 
do not decay at spatial inflnity. 

C. Resonant States 

We will now take a closer look at the energy interval < < /i^. Even though the 
reflection coefficient is one in this case (see Eq. (IT^ above), it is still possible to have 
resonances if Vcs{r*) has a "potential well" with a local minimum lying below /i^. In this 
case, plane waves emanating from r* — —00 may tunnel through the potential barrier into 
the potential well, and if they have the correct frequency their amplitude may build up 
through multiple reflections at the walls of the potential well. As we will show below, such 
resonances indeed occur for a discrete set of energy levels. 

Clearly, a necessary condition for this phenomena to take place is the existence of the 
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potential well. In order to find under which conditions such a potential well exists we first 
need to determine the critical points of the effective potential Vcs- It is in fact not difficult to 
show that such critical points correspond to the roots of the following third degree polynomial 
in r (notice that here we are working with r and not with r*): 

Mfi^r^ -£{i+l)r^ + 3M {f + i-l)r + 8M^ = . (14) 

This polynomial can be shown to have three real solutions if the following condition is 
verified: 

+ 7T^\/3(3£4 + Qp + 5£2 + 2£ + 3)3 , (15) 
zoo 

If M/i satisfies relation (ITSl) above, the potential V"cfr will have three critical points. One 
of them always corresponds to a negative value of r, so it is not of interest to our problem 
(negative values of r are unphysical). The other two critical points correspond to a minimum 
and a maximum of the effective potential. 

For the case /i = we recover the well-known solutions for which the maximum of the 
potential is located at: 



_^(. /. 14f 2 + 14^ + 9 1 

2 V 9^2(^+1)2 2£(£+l)'- ^ ' 

For £ — oo the maximum occurs at rmax = 3M, which represents the innermost stable 
circular orbit for massless classical particles. Any particle that starts with vanishing radial 
velocity in the region r < 3M will fall into the black hole. 

The region in parameter space where a potential well exists and resonances are possible 
is illustrated in Fig. |2]for different values of I (see also Fig. [T]). We will call this region the 
"resonance band", V^^™ < < min {VJ!^'^^, /i^}, although not all the states in that band 
will be resonant. Here V"g™™ and V^^^ are the local minimum and maximum, respectively, 
represented with dots in Fig. [H The cut-off in M/i that appears for each value of £ corre- 
sponds to a value of M/i large enough so that relation f|T5|) is no longer satisfied and the 
local minimum of the potential does not appear. Notice that such a cut-off increases with 
so that for a mass of the scalar field such that M/x < 1/4 the effective potential will have 
a minimum for all possible values of [f] We also want to stress the fact that we can have a 
^ Taking the mass of a typical black hole in the center of a galaxy to be of the order M ^ lO^Mo we obtain 
that the mass of the scalar field should be less than 10"^^ cV. Compare that result with the expected 



5o|, 



mass for the axion, 10 ^ to 10 eV 15(1, and with the mass of the SFDM models, 10 to 10 ^''^eV. 
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FIG. 2. Resonance band in parameter space. The cut-off in M ^ is reached for values so large that 
condition (jl5p is no longer satisfied. For 1 = this condition is M/i < 1/4, while for 1 = 1 one 



finds {Mfif < -9/32 + V20577/288 ^ 0.217 , which implies Mfi < 0.466. 



minimum in the effective potential for £ = 0, but not for /i = 0. This is the reason why the 
phenomenology of the problem we are analyzing here is quite different to that encountered 
in the usual Regge- Wheeler equation. 

Having found under which conditions there is a potential well, we now turn our attention 
to the existence of resonant modes. In order to have modes that are "almost trapped" in 
the potential well it is clear that one must ask for < < fi^. This will guarantee that 
the solution for the scalar field decays exponentially at spatial infinity and is "localized" 
close to the black hole, but of course the scalar field can still "escape" towards the black 
hole horizon. The situation we are considering is as follows: We are solving the stationary 
Klein-Gordon equation, so close to the horizon (where r* approaches — oo) we must have a 
combination of incoming and outgoing modes of the same amplitude. This implies that any 
scalar field that "escapes" from the potential well (via "quantum tunneling") is precisely 
compensated by the incoming field from the horizon. The situation is clearly not physical 
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since, first, we do not expect anything coming from the horizon, and second, the energy of 
such solutions is infinite as discussed above, but nevertheless these are the type of solutions 
we are considering here. 

Now, as already remarked, the frequency spectrum is continuous, so we will find a solution 
with the properties just described for any < < /i^. Clearly, if we want to find solutions 
that are in some sense "trapped" inside the potential well, we must also ask for uj^ to be in 
the resonance band, V^g™ < < min {V^^"", /i^}. The question now is what distinguishes 
the set of discrete resonant solutions from all other solutions in the continuous spectrum 
of the resonance band. There are several ways in which one can look for such resonant 
states, but here we have settled on a practical "intuitive" approach: Resonant states will 
correspond to those that "leak out" of the potential well very slowly, or in other words, they 
will correspond to those special solutions for which the ratio between the amplitude outside 
the potential well and the amplitude inside is a minimum. 

In order to have a clear picture of what we are talking about here, let us consider the 
particular case ^ = 1, M[i = 0.3, which clearly satisfies the conditions for a potential well to 
exist. Figure [U shows the effective potential M'^Vcs in that case. Notice that there is a clear 
minimum at r*/M ^ 25 for which M^V^™" ^ 0.085, while the maximum is at r*/M ^ 2.2 
with MVg^'^^ ^ 0.13. 

We can now solve the Schrodinger equation (Q numerically for an arbitrary value of 
< cj^ < yU^. The details for the numerical solver are not particularly important, it is 
sufficient to say that the solver starts at a large value of r*/M, for which Vcs ~ /i^, and 
where one can approximate the solution with a decaying exponential. One then solves the 
equation numerically from right to left. 

Figure [3] shows two solutions for particular values of Mu. The first plot corresponds 
to the solution for Mu = 0.295. Notice that for r*/M < the solution is oscillatory as 
expected, and that it penetrates only slightly into the potential well: The amplitude for 
r*/M < is much larger than the amplitude inside the well. This behavior is quite generic. 
The second plot, on the other hand, correspond to the solution for Mu = 0.29619. Notice 
that now the situation is quite the opposite: The amplitude inside the well is much larger 
than the amplitude for r*/M < (both solutions have been scaled so that the maximum 
overall amplitude is equal to 1). This means that the value Mu = 0.29619 is close to a 
resonant frequency; in fact it is close to the very first resonance, the one with no nodes 
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r*/M 

FIG. 3. Numerical solutions of the Schrodinger equation ^ for different values of Adu. The upper 
panel corresponds to Mu = 0.295, while the lower panel corresponds to Mlo = 0.29619. Both 
solutions have been scaled so that the maximum overall amplitude is equal to 1. 

inside the potential well. 

In order to locate the different resonant frequencies one can find the solution for all 
frequencies in the interval V^^"^ < cj^ < /i^, and simply plot the ratio between the amplitude 
Aout of the solution for r*/M < 0, and the amplitude Ai^ for the solution inside the potential 
well. Figure m is just such a plot for the case we are considering here. Notice that the ratio 
y4out/^in drops to vcry small values at a set of sharply defined discrete frequencies. Each 
subsequent resonant frequency can be easily shown (by just plotting them) to correspond 
to solutions with more nodes inside the potential well. As the value of u approaches the 
parameter /i the resonant frequencies become closer to each other, and in fact there seem 
to be an infinite number of them. Care must be taken when finding the numerical solution, 
since higher resolution is required as u approaches /i (both in r*/M and in u), while at the 
same time one needs to consider much larger values of r*/M since the solutions become so 
much wider (specific examples of such resonant solutions will be considered in the following 
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FIG. 4. Resonant frequencies. We show a log plot of the ratio ylout/^in between the amplitude 
Aout of the solution for r* /M < 0, and the amplitude A^^ for the solution inside the potential well 
r*/M > 10. 

sections) . 

III. LONG LASTING CONFIGURATIONS. NUMERICAL STUDY 

According to the arguments given in Section |IT1 above, we can not have stationary, lo- 
calized, physical scalar field configurations surrounding a Schwarzschild black hole: The 
stationary solutions, even the resonant ones, have a divergent energy due to their oscillatory 
behavior close to the horizon. However, one can always construct physical configura- 
tions that are arbitrarily close to the stationary solutions and that can survive the black 
hole for arbitrarily long times: Simply choose any (non-physical) stationary solution that 
decays at spatial infinity (i.e. < < yU^), and set it to zero by hand in the interval 
r G (2M, 2M + e), for some small parameter e > with dimensions of length. The resulting 

^ Moreover, such solutions are not compatible with the test field limit approximation. 
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configuration can be seen as a combination of tlie stationary solution plus a perturbation 
around the horizon, and it will be stationary for the domain of dependence of [2M + e, oo). 

For a given value of e, and some A > e, it will take a time of order t^{2M + A) = 
2Mln(A/e) + (A — e) for the perturbation to reach the radial coordinate r = 2M + A. That 
is, the standard fiat space-time result A — e, plus an extra term coming from the non-trivial 
space-time curvature, 2Mln(A/e). In the limit e — )■ that time can be made arbitrarily 
large. Since the energy density diverges as r — )■ 2M, one could naively think that we can not 
make e arbitrarily small if we want to remain within the test field approximation, but notice 
that, since the Klein-Gordon equation is linear in 0, for any given value of e one can always 
reduce the amplitude of the scalar field arbitrarily. But of course, that would yield scalar 
field distributions in which most of the energy is concentrated very close to the horizon. 

If we take 2M -|- A to be of the order of the size of a given scalar field distribution 
surrounding a Schwarzschild black hole, one would expect that t^{2M + A) could give us 
some insight into the lifetime of that particular configuration. f| It is worth noticing, however, 
that te works only as a lower bound, and some scalar field distributions could very well last 
much longer than that. In fact, we will show that initial data can be constructed in a 
particular way that produces long lasting, quasi-stationary scalar field distributions in a 
compact region outside the black hole. These initial data correspond to the resonant states, 
and we will devote most of the remaining of this paper to study them, although we will also 
consider a few other configurations for comparison. We will come back to the study of the 
characteristic time-scales associated to these configurations in Section IIII C 21 



A. Initial Data 



We will now construct pseudo-stationary configurations by truncating the stationary 

solutions as explained in the previous section. These configurations will then be used as 

initial data to be evolved numerically so that we can study the energy decay in a compact 

region outside the black hole. As already mentioned, the spectrum of stationary solutions 

is continuous, however we will show that pseudo-stationary initial data constructed from 

a particular discrete subset of stationary solutions, that of the resonant states, yields long 

^ In fact, for the configurations wc will be considering here (see Table this lifetime will not be particularly 

sensitive to the parameter e. If the width of the scalar field configuration is of order lO^A/, then one needs 

to take e < 10^^^ for the gravitational term 2A/ln(i?/e) to be relevant in the expression for te- 
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lasting configurations. 

As described above, we construct initial configurations that correspond to stationary 
solutions in the domain of dependence of the region r > 2M + e, with e < 0.05M . To do so, 
we integrate Eq. numerically in the region [R[^ = 2M+e, -Rout] using a shooting algorithm 
(see e.g. js^). As boundary condition we impose at r = i?out the relation u' = —\k\u with 
the constant k defined in Eq. (fT2l) to ensure that u{r) falls exponentially at large r (see 



also j46| for more details). For the left boundary condition we just impose u(-Rin) = 1, since 
the equation is homogeneous so u{r) can be rescaled afterwards. Finally, u{r) is simply set 
to for r < This introduces a discontinuity in the scalar field configuration, but the 
discontinuous jump is very small compared to the maximum absolute value of the scalar 
field and does not seem to introduce any numerical artifacts. 

Leaving the value of e aside, these initial configurations are characterized by the param- 
eters i and M/i, which parametrize the effective potential M^Vcs{r* / M) , together with Moj 
(see Eqs. ([9]) and (fTOj) above). Here we will restrict ourselves to the region < u;^ < /i^. 
The solutions with > yU^ are not localized close to the black hole and they will not be 
interesting for our purposes here. For the subspace with < oj^ < fj,"^ one can distinguish 
the following cases: 

1. Configurations constructed from stationary solutions in the resonance band, i.e. with 
V^f' <uj^ < min {V^^^^ /i^}. 

(a) Pseudo-stationary configurations constructed from resonant stationary solutions. 
We will call these pseudo-resonant states. 

(b) Pseudo-stationary configurations constructed from non-resonant stationary solu- 
tions. We will call these non-resonant states. 

2. Configurations constructed from stationary solutions outside the resonance band. 

From now on we will also refer to these configurations as type [Tal type [lb] and type |2l 
respectively. Here we are assuming that the effective potential does have a local minimum 
for some value of r*, i.e. Eq. ( IT5l) is fulfilled. If that is not the case, all the configurations 
will be of type 121 

Examples of type [la] and type [2] configurations are shown in Figure [5] where we plot 
the effective potential, together with two eigenvalues {Mu)^ (top panel), and the radial 
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energy density of the corresponding eigenfunctions (bottom panel). The Figure is shown 
in Eddington-Finkelstein coordinates, more appropriate for the numerical evolution (see 
Section nil Bp . For this particular case we have chosen i = 1 and Mfi = 0.3, so that Eq. (IT^ 
is fulfilled and a potential well is present (that case was also considered in Section [II CD . 
One would expect that such solutions will have an oscillatory behavior in regions where 
Vcs < uj"^, with most of the scalar field energy contained in those regions; and with a rapid 
decay outside. To help guide the eye, we added vertical lines to the plots at the locations 
where Ves = (classical turning points). We can see that most of the scalar field is located 
inside the classically permitted regions (Vcs < w^). In particular, the type [2] configuration, 
for which < V^g'"^, is located mostly to the left of the centrifugal barrier. A type [lb] 
configuration has not been shown here, but we can look at these configurations as somehow 
intermediate between type [la] and type [2] 

We will now concentrate on type [la] configurations which, as we will see in Section IIIIBl 
are the ones that show long lasting scalar field distributions at a compact region outside the 
black hole. We construct pseudo-resonant initial data configurations for different values of 
i and M/i, and for each pair of these parameters we study pseudo-resonant states up to the 
fifth mode. 

Table [T] summarizes some properties of the pseudo-resonant state initial data we have 
studied here. In what follows, -R99 is defined as the radius of a sphere containing 99% of the 
scalar field energy E at t = 0, and gives a rough idea of the "size" of each configuration. 
Note that -R99 increases as Mfi decreases and n increases. Some of these configurations are 
illustrated in Figures |6] and [7] Figure |6] shows the energy density distributions pe for the 
first five pseudo-resonant modes n = (1,2,3,4,5) for £ = 1 and Mfi = 0.2, while Figure [7] 
shows the first mode n = 1 pseudo-resonant configuration for different values of Mfi. 

Note that the eigenvalues w"^ of the resonant modes accumulate at w"^ = fi^ from below. 
It seems likely that there is in fact an infinite number of resonant modes n, although we have 

n 

not investigated in any detail if this is indeed the case. |j The possible existence of modes 
with arbitrary n would be interesting, since the energy distribution seems to spread further 
away from the black hole for larger tt,, implying that arbitrarily wide resonant modes might 
be constructed. 

^ Note that in the case of bound states in quantum mechanics, effective potentials with an asymptotic 
behavior for ?' — > 00 hke the one in Eq. (|10p usuahy have an infinite number of bound states. 
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FIG. 5. Top panel: Effective potential Vcg{r) for £ = I and Mfi = 0.3; and (MtOi)'^ for each of the 
solutions shown in the bottom panel. Bottom panel: Radial energy density Pe{i~) for the first (ni) 
pseudo-resonant mode (tvpe[Ta|) and a solution (labeled rix) with eigenvalue (labeled w^) below the 
local minimum of V^g (type[2|). An arbitrary rescaling was chosen for each in this plot in order 
to improve visualization. The vertical lines intersect the roots of V^s = w'^, a total of three for 
the pseudo-resonant mode (fii), and only one for the non-resonant mode (note that the root of the 
non-resonant mode is very close to the first root of the resonant mode, so the two corresponding 
vertical lines are indistinguishable in the figure). 

B. Numerical Evolution 

Having constructed a large set of initial data, we now study some properties of their 
evolution. The numerical evolution is performed in penetrating coordinates, hence covering 
a region that reaches inside the event horizon and avoiding the need to provide bound- 
ary conditions there. We use ingoing Eddington-Finkelstein coordinates {i,r,6,(p), with 
t := t — 2Mln(l — N[r)), and N{r) defined in ([2]). In these coordinates the spacetime line- 
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TABLE I. For the case i = 1, and different combination of parameters M^u and n we show: (i) The 
"width" of the initial energy distribution Rgg and (ii) The eigenvalue (Mcj)^. Notice that Rqq/M 
is given with three significant figure accuracy, while (Mw)^ is given with an accuracy up the last 
significant figure shown. All values shown correspond to pseudo-resonant (tvpe [Ta|) states. 
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element takes the form 



ds^ = -[l dP + dtdr + 1 + dr^ + r^dn^ . (17) 



This coordinate system covers a region that includes the black hole interior r G (0, oo), and is 
regular at the horizon. Notice that for the purposes of our simulations the Schwarzschild and 
Eddington-Finkelstein times can be considered equivalent since at a fixed radial coordinate 
r one can easily see that At = At. Because of this, and in order to simplify the notation, 
from now on we will drop the bar from the Eddington-Finkelstein time and will refer to it 
simply as t. 

In order to find numerical solutions to the Klein-Gordon equation, we define a set of first 
order derivatives, and obtain a system of evolution equations of the form dt"^ + Bdr'^ = S, 
with B a symmetric matrix and the vector formed from first derivatives in space and time 



of ipim- See references |46|, l52| for details. 

The evolution equations are solved numerically using second order finite differences in 
space, and evolving in time using a method of lines with a third order Runge-Kutta inte- 
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FIG. 6. Radial energy density pE at t = for different pseudo-resonant (type [Ta|) modes 
n = (1,2,3,4,5), for the case i = 1 and Mfj, = 0.2. Notice that the solution has been normal- 
ized so that we always have = 1. 



grator. Since we are using horizon penetrating coordinates we set the left boundary r^^i^ 
inside the event horizon, leaving the boundary conditions "free" at that point (by using 
one-sided spatial differences), since that region is causally disconnected from the exterior. 
On the other hand, the right boundary rmax is set far away and all incoming modes are 
set to zero there. The code passed the standard convergence and residual evaluation tests. 
Furthermore, two independent codes have been used in this work, and comparisons give 
identical results up to truncation error. 

The resolution used varies depending on the initial data. Typically, the region containing 
most of the scalar field (i.e. from to about -R99) is covered with at least 2,000, and often 
as much as 32,000, grid points. Then, the domain is chosen so that r^ax ~ 4i?99 or larger. 
Note that such high resolutions are needed in some cases in order to calculate the parameter 
s, defined below, with a relative error smaller than one. This is mainly due to its proximity 
to zero. 
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FIG. 7. Radial energy density pE for £ = 1 and the first pseudo-resonant (type [Ta|) mode n = 1, 
for configurations with different values of Mfi. The solution has been normalized so that = 1. 

C. Results 

We study the evolution of the scalar field ipim{t, r) for different initial data configurations. 
Although our attention is mainly focused on the evolution of the long lasting pseudo-resonant 
states, we also consider other configurations for comparison. We start by evaluating total 
energy loss and studying some spectral characteristics of the different configurations by 
means of a time Fourier analysis, and finish with more explicit considerations about how 
long can such configurations last. 

1. Energy Decay and Spectral Analysis 

In order to evaluate how long a scalar field configuration remains confined in a compact 
region we evaluate the conserved energy Ei^, defined in equation but we replace the 
upper limit of integration with a finite value R > Rgg, with typical values of R between 2Rgg 
and 10i?99. 

We begin by studying the evolution of the pseudo-resonant initial data (type [Taj) presented 
in table [B Figures [8] and [9] show the evolution of the integrated scalar field energy, E{t) (we 
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will drop the im sub- indexes from now on), for some of the parameters Mfi and Mun studied. 
Notice that the energy is plotted on a logarithmic scale, and shows a slow exponential decay 
of the form 

E{t) = Eo exp(-s t/M) , (18) 

with s constant, except for some very small oscillations that remain during the whole evolu- 
tion (these oscillations are noticeable in the overlay frames in the figures). Given the expo- 
nential decay that dominates the overall behavior, we can perform a linear fit of \n{E/Eo) 
as a function of t/M to calculate the parameter s. These values are shown in Table HT] below. 
In all the cases considered here no energy is seen to escape to infinity, but it rather falls 
into the black hole. In Section IIII C 21 below we will evaluate explicitly the characteristic 
times for these configurations. The results shown on the table correspond to the evolution of 
pseudo-stationary initial data constructed from the discrete spectrum of resonant stationary 
states (type [Ta|) . 
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FIG. 8. Energy of the scalar field vs. time for the evolution of pseudo-resonant (type [Ta|) initial 
data with i = 1, M/i = 0.35 and modes n = (1, 2, 3, 4, 5). As can be clearly seen some evolutions 
last longer than others, this is because we chose to run each case for a time length approximately 
proportional to the size of the scalar field distribution, Rgg. 

We have also evolved a couple of type [lb] configurations, that is, non-resonant pseudo- 
stationary states. Such non-resonant states do not evolve in a quasi-stationary fashion, 
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FIG. 9. Energy of the scalar field vs. time for the evolution of the initial data corresponding to 
the first (n = 1) pseudo-resonant (type [Taj) mode with i = 1, and M/i = 0.2, 0.25, 0.3 and 0.35. 



TABLE II. For the case with i = 1, and each combination of parameters Mfi and n, we show the 
slope s of a linear fit of ln{E{t)/EQ), as described in the text. The error in the slope is estimated 
either as one half the difference of the values obtained with the two higher resolutions used in each 
case, or as the error obtained from the linear regression, whichever is larger. 
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instead they very quickly lose a significant portion of their energy, some falls into the black 
hole while some is radiated toward infinity. It is only after quite some time has passed that 
they reach a quasi-stationary state, and with much less than their original energy. Once 
they are in the quasi- stationary regime, the energy falloff is very similar to that of the 
pseudo-resonant states. Each non-resonant state seems to evolve, at least after some time, 
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ClS Sb combination of pseudo-resonant modes. 

To understand this behavior better we perform a spectral analysis. We calculate the 
discrete Fourier transform in time of the numerical scalar field evolutions. More precisely, 
we calculate the magnitude of the discrete Fourier transform in time of the field ip{t,r) 
at a fixed spatial point (or points) r = rj, located approximately at a local extremum of 
ip{t = 0,r), and pay attention to possible differences derived from using different sample 
points Vj. The explicit expression is 



where A is a normalization constant and tp are the discrete time values. 

The Fourier transform of a (type [Ta|) pseudo-resonant state is shown in Figure [TUl and 
that of two (type [Tb|) non-resonant ones is shown in Figures [TT] and [T2i For this example 
we have chosen to show an n = 2 pseudo-resonant state and two non-resonant states with 
frequencies Ux between the first and second pseudo-resonant modes (wi < Ux < u}2)- These 
two frequencies have been denoted with an x in the overlay figures, while numbers from 1 
to 5 denote the first five resonant frequencies. All cases correspond to £ = 1 and Mfi = 0.3. 

Figures [TUl and ITT] show the Fourier transform at a single sample point ri. While the 
results shown in these two figures are almost independent of ri, for the case considered in 
Figure [12] there were some small differences depending on whether the sample point rj was 
chosen close to the first or the second local extremum, so we show both results in that case. 
Note however that the two curves in Figure [12] are qualitatively identical, and the peaks are 
located at exactly the same frequencies, the conclusions we will make from this figure are 
then independent of the sample point chosen. 

Having mentioned a few technical details about the spectral analysis, we now make some 
observations. Consider first the spectrum of the pseudo-resonant mode with frequency U2 
shown in Figure [10] As might be expected for this configuration, we see a single and well 
defined peak at exactly uj2- Note however that this result is not so obvious a priori, since 
we are actually evolving a truncated version of the actual stationary state. Consider next 
the non-resonant configurations. Figures [TT] and [T^ with frequencies u^- We note that no 
peaks are present in the spectrum at the frequency Ux, instead one can see clear peaks at 
different resonant frequencies Un- Note also that, although not strictly true for all cases, 
the peaks tend to be higher the closer they are to Ux- It would seem that, after some 




(19) 
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FIG. 10. Discrete Fourier transform in time vs. frequency for the evolution of type [la] data with 
M/i = 0.3, i = 1, and frequency W2 (corresponding to the second pseudo-resonant mode). The 
overlay figure shows in more detail the region contained in the dashed line rectangle. Numbers 
going from 1 to 5 in the overlay figure denote the first five resonant frequencies. 

transient, these non- resonant states evolve as a combination of pseudo- resonant states. In 
that sense, it seems that the pseudo-resonant states could play a similar role to that played 
by the quasi-normal modes in the study of perturbed black holes. This point should be 
investigated further, but it is outside the scope of the present paper. 

The Fourier transforms were calculated integrating over a time interval that goes from 
t = to t ~ lO^M. Such a large interval is necessary to ensure that the final result 
has enough frequency resolution to distinguish between the different resonant frequencies 
Un, which are very close to each other. This fact also implies that we cannot integrate 
over shorter intervals in order to evaluate possible time variations of the spectrum at the 
beginning of the evolution, before reaching a quasi- stationary state. 

2. Characteristic Timescales 

We now turn our attention back to the evolution of pseudo-resonant (type [Taj) configura- 
tions. Based on the results outlined in table [Til we can be more explicit about how long this 
states can last. If the exponential regime observed is sustained for all times, a characteristic 
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FIG. 11. Discrete Fourier transform in time for the evolution of type flbl data with Mfi = 0.3, i = 1, 
and frequency Wx = 0.29664794. The overlay figure shows in more detail the region contained in 
the dashed line rectangle. Numbers going from 1 to 5 in the overlay figure denote the first five 
resonant frequencies, while the x denotes Wx- 

time of these configurations, their half-hfe time, will be given by ti/2 := ln(2)M/s, with s the 
parameter for the exponential decay of Eq. f ll8l) . Note that in some of the cases considered 
we were able to evolve far beyond ti/2 thus corroborating, at least for those parameters, that 
the exponential regime is indeed sustained (see Fig. [8]). 

Changing back from geometric units, and considering a black hole mass of IO^Mq, we 
get a characteristic time of order 

680 

ti/2 ~ seconds . (20) 

s 

Replacing, for instance, the value s = 7 x 10~^ that corresponds to ^ = 1, M/i = 0.2, n = 4 
(see Table [IT|) . we obtain ti/2 ~ 3, 000 years. This is still a very small time when compared to 
cosmological time-scales, however, values of the parameter /i motivated by dark matter scalar 
field models correspond to hfi ~ 10~^^eV in physical units (see Introduction), or equivalently 
Mfi ~ 10~^ in geometric units, and noticing how fast the parameter s seems to decrease 
with decreasing M/x (see Table HTl) . one might expect that configurations with values of M/j, 
motivated by dark matter scalar field models will have half-life times of cosmological scales. 
We were, however, unable to study cases with such small so the preceding argument is 
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FIG. 12. Discrete Fourier transform in time for the evolution of type flbl data with M/x = 0.3, i = 1, 
and frequency Wx = 0.29782545. The overlay figure shows in more detail the region contained in 
the dashed line rectangle. Numbers going from 1 to 5 denote the first five resonant frequencies, 
respectively, while the x denotes Wy^. The solid line curve (sample 1) corresponds to a sample point 
located close to the first extremum, while the dashed line curve (sample 2) correspond to one close 
to the second extremum. 

purely hypothetical. 

Note that extrapolating to values of Mfi of order 10~^ from the few values given in 
Table HTl which are of order 10~^, would be extremely inaccurate. If, however, we were to 
proceed anyway with such an extrapolation, we would obtain s of about 6 x lO^^'^, implying 
a value of ti/2 of about 10^ years. Compare this value with the expected age of the universe, 
10^° years. U 

Finally, for the purpose of comparison with the long lasting f|Ta|) pseudo-resonant con- 
figurations, we will consider the evolution of configurations of type ([2]). In the solid-line 
curve of Fig. [13] we show the half-life time ti/2 for both pseudo-resonant states and states 
outside the resonance band for the case £ = 1 and different values of M/i. A change in the 
slope can be seen as soon as the critical mass Mfi = 0.466 is reached, showing that ti/2 
is indeed larger for type f|Ta|l states than for type ([2]) states. Here ti/2 is expressed in the 
^ In the parameter range studied, s decreases slightly faster than an exponential with decreasing fi. For the 
extrapolation we assume an exponential that fits the points with smaller M/i, Alfi = 0.2 and Mfi ~ 0.35. 



27 



10 



a 10 

H 



10^ 



0.1 



I ' ' ' ' I ' ' ' ' I 

Resonant region 



Non-resonant region 




0.2 



0.3 



0.4 



0.5 



0.6 



10 



10 



10^ 8 



a 

H 



10 



10 



FIG. 13. Solid line: Half-life time of different scalar field configurations, ti/2 ■= ln(2)M/s. Dashed 
line: "Collapse time" , defined as the time in which a pressure-less gas sphere of size Rgg /2 collapses 



in the gravitational field of the black hole, tc := a/ Rgg/8GM. Dotted line: Characteristic time 
defined at the beginning of Section Hill te{2M + Rgg/l). The left Y-axis shows these characteristic 
times in units of M, while the right Y-axis is expressed in seconds (assuming M = 10^ Mq). 



left vertical axis in units of M, while on the right vertical axis it is expressed in seconds 
(assuming M = 10^ Mq). 

A quick inspection of Table [I] shows that pseudo-resonant states have a Rgg larger than 
that of non-resonant states, i.e. they are wider. One can argue that this is the reason why 
pseudo-resonant states last longer. In order to show that this is not the case, we can make a 
comparison with some other characteristic times. For instance, we can compute the collapse 
time tc, which is the time in which a pressure-less gas configuration of size a collapses by 
its own gravity. This is defined as tc := a/ct'^/GM^, where Mg is the mass of the gas. In 
our test field approximation, Mg will be the mass of the black hole Mg = M. In order to 
compare tc with ti/2 we found appropriate to choose a as Rgg/2. This time is shown in 
the dashed-line curve of Fig. [131 As can be clearly seen in this figure, tc is always smaller 
than ti/2, and for the particular case of pseudo- resonant states ti/2 is orders of magnitude 
larger than tc- Finally, for comparison, in the Figure we also show (dotted-line curve) the 
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characteristic time defined at the beginning of Section lllli 



IV. CONCLUSIONS 

We investigated the existence of long-lasting scalar field configurations surrounding a 
black hole. A good motivation for this study is the possibility that super-massive black 
holes at galactic centers may represent a serious threat to the scalar field dark matter 
models. This is because such black holes may in principle swallow an hypothetical scalar 
field halo in cosmologically short times, which is what in general will happen with arbitrary 
scalar field distributions. However, we were able to find particularly long-lasting scalar field 
configurations around a black hole. 

As a first step, we have considered a relatively simple model and looked for solutions to 
the Klein-Gordon equation on a Schwarzschild space-time background. Although we do find 
stationary solutions for the scalar field, we show that they are unphysical, in the sense that 
their energy density integrates to infinity in a compact region just outside the event horizon. 
However, we were able to find long-lasting, quasi-stationary solutions of finite energy. This is 
done by evolving initial data that was constructed by modifying slightly a particular subset 
of the (unphysical) stationary solutions. The solutions found show as an overall behavior 
an exponential energy decay, caused by scalar field leaking into the black hole, that in some 
cases can be very slow. 

The stationary solutions are obtained by solving a time-independent Schrodinger-like 
equation with an effective potential. Hence, they can be characterized by the properties of 
that potential. This fact is strictly true for the stationary solutions, but interestingly we 
find that the (physical) quasi-stationary solutions, for which the Schrodinger-like equation 
no longer holds, are also characterized by properties of that same effective potential. The 
cases of interest are those in which the effective potential has a local minimum (a potential 
well). The existence of this minimum depends solely on the combination of parameters M^x 
and £. Although none of the possible forms of the effective potential allow for bound states, 
the existence of the mentioned well is enough to allow for the so-called resonant states. 
These states are the ones that have proven useful for constructing initial data that give rise 
to long-lasting quasi-stationary configurations of finite energy. 

It may be objected that in order to obtain the mentioned quasi-stationary solutions 
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one has to construct very particular initial data, a situation that might be very unlikely 
in nature. However, it seems that the crucial factor is only the existence of the potential 
well. We showed that even when starting with modified stationary solutions that are not 
resonant, after an initial abrupt energy loss, the late time behavior observed is very similar 
to that of the resonant quasi- stationary solutions. In fact, these solutions seem to evolve 
as a combination of the resonant modes. Note that the value of for scalar field dark 
matter models is expected to be given approximately by = 10~^^eV in physical units, or 
M/i ~ 10"^ in geometric units for a IO^Mq black hole, which gives rise to effective potentials 
with a local minimum for all values of i. 

When evaluating the characteristic times of our solutions we find that the longest lasting 
configurations last for times of the order of thousands of years. Although this is cosmologi- 
cally a very short time we must note that, for technical reasons, we were only able to study 
cases with relatively large values of Mfi, of order 10~^. Noting how fast the characteris- 
tic times seem to increase with decreasing fi, it would not be unreasonable to expect that 
configurations with M/i ~ 10~^ could last for cosmological time-scales. One could argue in 
the same lines about the "size" of the scalar field distributions: The configurations studied 
extend only up to a radius Rgg ~ 2, DOOM, but again, this number increases rapidly with 
decreasing /x. 

A few aspects of our study are still open for some improvement, and will be addressed in 
future works. First, a self-gravitating scalar field may be considered. Due to the extremely 
diluted nature of dark matter halos, the test field (Cowling) approximation is very accurate 
to some degree. However, given the halos large spatial extension, one expects that important 
differences may appear globally as a cumulative effect. Second, much smaller values of the 
parameter /i, and much larger scalar field distributions, will be needed to do a more realistic 
representation of dark matter halos. The main difficulty in dealing with such configurations 
is handling numerically the very different scales, going from the regions close to the black 
hole to galactic-size scales. Mesh refinement could be used to solve this difficulty. Third, 
besides studying possible quasi-stationary or long-lasting configurations with an already 
existing black hole, it would be interesting to consider more dynamical scenarios such as 
the formation and/or growth of the black hole and the possibility of survival of the scalar 
field afterwards. We plan to address these and other issues in future works. However, the 
results presented here already seem to indicate that it is indeed possible for scalar field halos 
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around super-massive black holes to survive for cosmological time-scales. 
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